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Abstract. We prove short-time existence for the Einstein-Euler-Entropy 
system for non-isentropic fluids with data in uniformly local Sobolev 
spaces. The cases of compact as well as non-compact Cauchy surfaces are 
covered. The method employed uses a Lagrangian description of the fluid 
flow which is based on techniques developed by Friedrich, hence provid- 
ing a completely different proof of earlier results of Choquet-Bruhat 
and Lichnerowicz. This new proof is specially suited for applications to 
self-gravitating fluid bodies. 
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1. Introduction. 

The Einstein equations have been a source of several interesting problems 
in Physics, Analysis and Geometry. Despite the great deal of work that has 
been devoted to them, with many success stories, several important questions 
remain (see [13, 17, 50] for an account of what is currently known and some 
directions of future research) . One of them is finding a satisfactory theory of 
isolated systems, such as stars, both from a perspective of the time develop- 
ment of the space-time, as well as from the point of view of the geometry 
induced on a space-like three surface. To quote Rendall, "of the physical situ- 
ations which can be described by the general theory of relativity, those which 
are at the present most accessible to observation are the isolated systems. In 
fact, all existing tests of Einstein's equations concern such situations. It is 
therefore important to have a theory of these systems which is as complete 
as possible, not only in terms of the range of phenomena which are covered 
but also with respect to logical and mathematical solidity" [45]. 

Stars are the prototypes of isolated systems. They are typically modeled 
by considering a region of space filled with a fluid and separated from an 
exterior that corresponds to vacuum. The properties of the fluid, such as 
perfection, viscosity, charge, etc., depend on the particular situation one is 
interested in. The dynamics of the fluid region is then described by Einstein 
equations coupled to matter, whereas vacuum Einstein equations hold on the 
complement of this set. From the point of view of the Cauchy problem, which 
will be the main case of interest in this work, there are two primary questions 
to be addressed. 

(i) First, the solvability of Einstein equations in the two different sit- 
uations of interest, i.e., coupled to fluid sources and in vacuum, should be 
addressed. Since the short-time existence for vacuum Einstein equations is 
well understood (see e.g. [30]), this leaves us with the coupling to matter. A 
specific matter model has then to be chosen, and one of the most customary 
choices, sufficient for many applications, is that of a perfect fluid [2], in which 
case Einstein equations are coupled to the (relativistic) Euler equations. The 
well-posedncss 1 of this system was proven by Choquet-Bruhat [10] and ex- 
tended by Lichnerowicz [34] to include entropy — in which case the resulting 
system will be called the Einstein-Eulcr-Entropy system, whose equations are 
stated in section 2.2. 

(ii) The second question is more delicate and consists of trying to bring 
together the two different scenarios described in the previous paragraph, 
namely, vacuum and coupling to matter. More precisely, we attempt to for- 
mulate and solve the Cauchy problem with an interior region — thought of 
as the star — governed by the Einstein-Euler (or Einstein-Euler-Entropy) 
equations and an exterior one that evolves according to the vacuum Einstein 
equations. This corresponds to a genuine free-boundary problem in that the 
boundary of the star cannot be prescribed for time t > 0, being rather a 



1 "Well-posedness" should be understood here in the context of General Relativity, where 
uniqueness is meant in a geometric sense, i.e., up to isometries. 
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dynamic quantity that has to be determined from the evolution equations. 
The star boundary at time zero is obtained as the boundary of the support 
f2 of the initial matter density go, which is a function on the initial Cauchy 
surface typically of the form 



Solving the Cauchy problem then requires a refined analysis of the bound- 
ary behavior of the quantities involved, where the change from G a p = T a p 
to G a p = causes severe technical difficulties (here G a p and T Q( g are the 
Einstein and stress-energy tensor, respectively). 

Lindbom [38] has proven that any isolated static fluid stellar model 
ought to be spherically symmetric, generalizing a classical result of Carleman 
and Lichtenstein for Newtonian fluids [7, 36, 37]. More precisely, he has shown 
that a static asymptotically flat space-time that contains only a uniform 
density perfect fluid confined to a spatially compact world tube is necessarily 
spherical symmetric. 

Under the assumption of spherical symmetry, it is possible to deal with 
many of the technicalities that arise, and a number of satisfactory and gen- 
eral results have been obtained. Rendall and Schmidt have proven existence 
and uniqueness results for global solutions of the Einstein-Eulcr system [47] . 
They also gave a detailed account of how the properties of these solutions 
depend on features of the equation of state. Kind and Ehlcrs have treated the 
mixed initial-boundary value problem in [33] , where they have also given nec- 
essary and sufficient conditions for attaching the solutions they constructed 
to a Schwarzschild space-time. Makino [41] refined the results of Rendall and 
Schmidt by providing a general criterion for the equation of state, which en- 
sures that the model has finite radius (and therefore finite mass) . This result 
is used to study the linear stability of the equations of motion. Regarding the 
relation between mass and radius, in [31] the authors derive, among other 
results, interesting mass-radius theorems. A very extensive treatment of the 
Cauchy problem for spherically symmetric data is given in the work of Groah, 
Smoller, and Temple [29]. These works do not exhaust all the results known 
in spherical symmetry; see the references in the above papers for more details. 

If we drop the hypothesis of spherical symmetry, however, not much 
is known about the well-posedness of the Cauchy problem when the matter 
density is allowed to vanish outside a compact set, as in (1.1) — a situation 
generally referred to as a "fluid body". One could attempt to change go to 
a function that is always positive, but decays sufficiently fast and is very 
small in the region outside Q. Unfortunately, this does not improve things 
considerably because generally the time-span of solutions cannot be shown to 
be uniform, and as a result, the domain of definition of the solution contains 
no space-like slice t = constant (except obviously the t — slice). Another 
problem which arises in this context is that the usual conformal method for 
solving the constraint equations [35, 62, 14] cannot be applied. 




(1.1) 
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Following the ideas of Makino [40] , who studied non-relativistic 2 gaseous 
stars, Rendall has shown that the above problems can be circumvented for 
fluids obeying certain equations of state [46] — although, as the author him- 
self points out, the solutions obtained in these cases have the undesirable 
property of not including statically spherically symmetric space-times, and 
the restrictions on the equation of state are too strong. Initial data sets for the 
Einstcin-Euler equations describing a fluid body have been successfully con- 
structed by Brauer and Karp [3] and Dain and Nagy [16], but well-posedncss 
of the system with such prescribed initial values has not yet been demon- 
strated. Therefore, the solution to the Cauchy problem for a self-gravitating 
isolated fluid body is still largely open (although some recent results of Brauer 
and Karp [4, 5, 6] address directly some of the technical issues we referred 
to). 

An important step forward has been achieved by Friedrich in [21]. Using 
a Lagrangian description of the fluid motion, he has been able to derive 
a set of reduced equations which form a first order symmetric hyperbolic 
system. While it has been known that the vacuum Einstein equations can 
be cast in such a form since the work of Fischer and Marsden [20], and 
the Einstein-Euler system had also been investigated within the formalism 
of first order symmetric hyperbolic systems in the aforementioned work of 
Rendall, what makes Friedrich's construction particularly attractive is the use 
of the Lagrangian description of the fluid flow, as it is known that Lagrangian 
coordinates are uniquely suited for treating free boundary problems involving 
the non-relativistic Euler equations, and have in fact been employed to a great 
success to study them [18, 15, 55]. 

It remains to be seen whether Friedrich's ideas will lead to satisfactory 
existence theorems for the Cauchy problem for fluid bodies. But for such a 
project to be successful, one needs to be able to carry out step (i) as explained 
above, namely, to use the reduced system of [21] to separately solve Einstein 
equations in the cases of vacuum and coupled to fluid matter. Only then the 
passage across the boundary from the Einstcin-Eulcr-Entropy to the vacuum 
equations can be analyzed. Notice that in the case of vacuum, the notion of 
a Lagrangian description is somewhat artificial, but it can be given meaning 
by the choice of a time-like vector field which should agree with the fluid 
four-velocity once matter is introduced. 

Well-posedness of the vacuum Einstein equations using the formalism 
of [21] has been established in [22, 23] (see also [24. 25, 26]). This leaves us 
with the question of well-posedncss of the Einstcin-Eulcr-Entropy system. In 
other words, we seek a solution to the following: 

Establish a well-posedness result for the Einstein- Euler- Entropy system 
with a general equation of state, by writing the set of equations as a first order 
symmetric hyperbolic system via a Lagrangian description of the fluid flow, as 



2 It should be noticed that even self- gravitating Newtonian fluid bodies are not well un- 
derstood. In fact, with the exception of the viscous case studied by Secchi [51, 52, 53], not 
many results besides that of Makino [40] seem to be available. 



Einstein-Euler-Entropy system 



■3 



in [21], and under the assumption that the fluid fills the entire initial Cauchy 
surface S, i.e. Qq > c > 0, with possibly additional hypotheses consistent with 
physical requirements. 

This is the problem addressed and solved in this work. Some terminology 
is needed before we can make a precise statement — see theorem 2.6 in 
section 2.2. We stress that theorem 2.6 had been proven much earlier by 
Choquct-Bruhat [10] and Lichnerowicz [34]. What is new is the method of 
proof, employing the Lagrangian description of the fluid flow. It should also 
be noticed that, while the equations we use are essentially those of [21] and 
[28], to the best of our knowledge, a complete proof of well-posedness relying 
on Friedrich's techniques is not available in the literature. 

1.1. Outline of the paper and notation. 

The paper is organized as follows. In section 2, we recall several basic defi- 
nitions, fix our notation, introduce the main equations and hypotheses, and 
state the main theorem. In section 3, we introduce the frame formalism and 
gauge conditions that constitute the basis of Friedrich's method. The gauged 
or reduced Einstein-Euler-Entropy system is also presented in this section. 
The proof of the main theorem is carried out in section 4. It consists of three 
parts: determination of the initial data (section 4.1), well-posedness of the 
reduced system (section 4.2), and the "propagation of the gauge" (section 
4.3). This last step is what guarantees that a solution of the Einstein-Euler- 
Entropy system in a particular gauge — the reduced system — yields a 
solution to the original set of equations. Finally, in section 5, we make some 
closing remarks. 

We have chosen to present our arguments in a logical rather than con- 
structive order. This means that instead of starting by showing how the 
reduced equations are obtained from the Einstein-Euler-Entropy system by 
means of a gauge choice, we first state the reduced equations, then derive its 
solutions, and finally show that they correspond to solutions of the original 
equations of motion. This particular order is adopted because the equations 
we use in both the reduced and the propagation of the gauge systems are 
equivalent to those of [21] and [28], although we shall write them in a slightly 
different fashion. Hence, as they have been considered before, there is no 
immediate need to show how the reduced equations are extracted from the 
original ones by a suitable choice of gauge. This procedure is briefly presented 
nonetheless in the appendix for the reader's convenience. Those not familiar 
with the work [21] are encouraged to read appendix B prior to section 4. 

Notation 1.1. The Sobolev space between manifolds M and N will be de- 
noted by H S (M, N). When no confusion can arise, the reference to N will 
be suppressed, and H S (M,N) will be written H S (M), or even H s when M 
is clear from the context. denotes the uniformly local Sobolev spaces, 
whose definition is recalled in appendix A. 
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2. The basic setting and the main result. 

Our main object of study is a four-dimensional Lorcntzian manifold (A4,g), 
called a space-time. In General Relativity, we are interested in the sub-class 
of Lorentzian manifolds where Einstein equations are satisfied. We recall that 
these are 

Ra/3 - 2 R 9a/3 = ^T a/3 , (2.1) 

where i? Q/ g and R are respectively the Ricci and scalar curvature of the metric 
g, T a f3 is the stress-energy tensor which encodes information about the matter 
fields, and ,J/if = 87r§, with c being the speed of light and G Newton's 
constant 3 . A space-time is called an Einsteinian space-time when (2.1) is 
satisfied. M. will always be assumed to be oriented and time-oriented. 

The left-hand side of (2.1) is divergence free as a consequence of the 
Bianchi identities. Hence, regardless of the particular matter model which is 
considered, the stress energy tensor has to satisfy 

V a T aP = 0, (2.2) 

which then gives equations of motion for the matter fields (see section 2.2 
below); V in (2.2) denotes the Levi-Civita connection associated with g. 
Equation (2.2) is sometimes referred to as the local law of momentum and 
energy conservation. 

Convention 2.1. Throughout the paper we shall adopt the convention (+ — 
) for the metric. 

Definition 2.2. A fluid source is a triple (U, g, u), where (U, g) is a domain of 
a space-time and u a time-like vector field on U of unit norm. Physically, the 
trajectories of u represent the flow lines of matter. Given a fluid source, the 
stress- energy tensor of a perfect fluid is given by 

T a i3 = (p + g)u a up - pg a ji, (2.3) 

where p, called the pressure of the fluid, and g, called the density of the fluid, 
are non-negative real valued functions. A perfect fluid source is a fluid source 
together with a stress-energy tensor given by (2.3). 

Perfect fluid sources, or perfect fluids for short, are often used to study 
space-times where a continuous distribution of matter exists (see e.g. [60]). 
The assumption that matter is described by a stress-energy tensor of the 
form (2.3) means that no dissipation of any sort is present; in particular one 
neglects possible effects due to heat conduction, viscosity or shear stresses 4 . 



3 As it is customary to drop the factor 87r, and we adopt units where G = c = 1, we could 
set = 1. But since we will be working with a lesser familiar formalism, it is useful to 
keep the constant Jf , in which case it is always possible to set J€ = in order to compare 
the resulting equations with the well studied vacuum case. 

4 In fact, the (here postulated) stress-energy tensor (2.3), can be motivated from basic 
physical principles, which can be taken as characterizing ideal fluids; see e.g. [2, 19]. 
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For perfect fluids, equation (2.2) becomes 




(2.4) 



+ (P + Q)u a ^aUp + upu a \7 a p - \7pp = 0. 

Taking the inner product of (2.4) with u and using u a u a = 1 (which also 
implies u"V pu a = 0), gives 



known as the conservation of energy or continuity equation. Then, using (2.5) 
into (2.4) produces 



known as conservation of momentum equation. Together, equations (2.5) and 
(2.6) are known as the Euler equations for a relativistic fluid. 

In physically relevant models, the functions p and g are not independent 
but are related by what is called an equation of state, which is an (usually 
smooth) invcrtible function 



Fluids where (2.7) is satisfied arc called barotropic fluids, with the particular 
case p = called dust or pressure-free matter. Physical situations of interest 
where barotropic fluids are employed include some models of cold (more pre- 
cisely, zero temperature) matter, such as completely degenerate cold neutron 
gases (which are used to model nuclear matter in the interior of neutron stars) 
[8, 2, 54]; the so-called ultra-relativistic fluids, i.e., fluids in thermal equilib- 
rium where the energy (which in relativistic terms is described by g due to 
the equivalence of mass and energy) is largely dominated by radiation [2, 60]; 
and fluids of electron-positron pairs [60] . Barotropic fluids are also important 
in Cosmology, where several models of the early universe assume that the 
distribution of matter is described by an ultra-relativistic fluid [8, 2, 60, 61]. 

As usual in General Relativity, space-time itself is a dynamic quantity, 
therefore its geometry, encoded in the metric g, as well as the the dynamics 
of other fields present on space-time, have to be determined as solutions to 
the Einstein equations (2.1) coupled to matter via (2.2). In particular, the 
fluid source (U, g, u) is not a given, but has to arise from the solutions of the 
coupled system. The Einstein- Euler system for barotropic fluids, or barotropic 
Einstein- Euler system, is given by equations (2.1), (2.5), (2.6), and (2.7), with 
T a p given by (2.3), and subject to the constraint u a u a = 1. The unknowns to 
be determined are the metric g, the four-velocity u, and the matter density 

The Einstcin-Eulcr system for barotropic fluids has been studied by 
many authors. In fact, most of the results cited in the introduction deal with 
this case. The interested reader can consult the papers [16, 3, 46, 31, 29, 41, 
47] and the monographes [34, 2, 8], as well as the references therein. 

Similar to the well-known case of vacuum Einstein equations [9], when 
addressing the solvability of the barotropic Einstein-Euler system, one has 



(p + g)V a u a +u a V a g = 0, 



(2.5) 



(p + g)u a V a u p + U(su a V a p - Vpp = 0, 



(2.6) 



P = p{g). 



(2.7) 
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to investigate its constraints and make suitable choices for the spaces where 
solutions will be sought. Since we shall deal with the more general case of 
fluids that are not necessarily barotropic, we postpone this discussion for the 
time being, turning our attention first to some thermodynamic tools that will 
be necessary in the sequel. We shall return to the barotropic case in section 
5. 

2.1. Thermodynamic properties of perfect fluids. 

There are important scenarios where it is known that the pressure is not de- 
termined by the matter density only. These include the so-called polytropic 
fluids, which are used in several stellar models (see e.g. [8, 2, 60, 63] and 
references therein). In such cases, (2.7) has to be replaced by a more gen- 
eral equation of state involving other physical quantities, which are usually 
assumed to be of thermodynamic nature, as we now describe. 

As in many situations in General Relativity, it is important to identify 
those quantities measured by a local inertial observer. Let r be the rest mass 
(or energy) 5 density 6 , defined as the mass density measured in the local rest 
frame. It is assumed that this quantity obeys the conservation law 

y a (ru a ) = 0, (2.8) 

which states that mass is (locally) conserved. Notice that on dimensional 
grounds, - is the specific (i.e., per unit of mass) volume. The difference 
between the mass density g and the rest mass density r is by definition the 
specific internal energy e (as measured in the local rest frame): 

e = r(l + e). (2.9) 

Therefore, as p, r and g, e is a real valued function on U C J\A. In the context 
of rclativistic fluids, the relation (2.9) has been first introduced by Taub [57] 
and has been widely used since [28, 8, 2, 34]. 

For non-barotropic fluids, further relations among the variables of the 
problem have to be introduced in order to have a well-determined system of 
equations. It is natural to assume that the first law of thermodynamics holds, 
i.e., 

de = K ds — pdv, 

where K is the absolute temperature, s the specific entropy, and v the specific 
volume. These are all no n- negative real valued functions onM. In light of (2.9) 
and v = the first law can be written as 7 

dg = ?-t-idr + rK ds. (2.10) 



5 Due to the equivalence of mass and energy, and our choice of units with c = 1, we shall 
use the terms mass and energy interchangeably. 

6 Also called particle number density [28], baryon number density [2] or yet proper material 
density [34]. 

7 The case of interest in this paper is when the matter density does not vanish, hence - is 
well defined because of (2.9). The case of vanishing q is, however, important, as explained 
in the introduction. 
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We now turn to the generalization of (2.7), i.e., to the appropriate equation 
of state that has to be provided. Although we presently have seven thermo- 
dynamic variables, namely, p, g, r, e, K, v and s, the above equations imply 
relations among them. In fact, for a perfect fluid, only two of such quan- 
tities are independent [2], with the remaining ones determined by relations 
depending exclusively on the nature of the fluid. On physical grounds, we 
should assume that such relations are invertible, what renders the question 
of which two thcrmodynamical quantities are the independent ones a matter 
of preference. We shall assume henceforth that r and s are independent. We 
thus postulate an equation of state of the form 

g=&(r,a), (2.11) 

where is a given smooth function, invertible in the sense that we can solve 
for r = r(g, s) and s = g{g,r). From (2.10) and (2.11) we then obtain 

P = r%~Q, (2-12) 



and 



K=±* (2.13) 
r os 

The validity of (2.10) along with (2.5) implies 

Ku a V a s = -\(p + g)V a {ru a ). (2.14) 

Therefore, in light of (2.8), and assuming that the temperature is not zero 
(which is consistent with (2.13) and (2.11)), 

u a y a s = 0. (2.15) 

In other words, the entropy is conserved along the flow lines of the fluid. 
When (2.15) holds, the motion of the fluid is said to be locally adiabatic. A 
fluid is said to be isentropic if s =constant, and non-sentropic otherwise. 

Remark 2.3. If (2.8) is not assumed, then from (2.14) it only follows that 

V Q (ru Q ) <0, 

and 

u a y aS > o. 

Although these two inequalities have a clear physical interpretation — from 
the point of view of inertial observers the rest mass cannot increase and the 
entropy cannot decrease — if (2.8), and therefore (2.15), is not assumed, the 
motion of the fluid is underdetermined. 

The first law of thermodynamics, equation (2.10), will always be as- 
sumed to hold, therefore, because of (2.14), we shall work with (2.15) rather 
than (2.8) to form our system of equations. 
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2.2. The Einstein-Euler-Entropy system. 

The Einstein-Euler-Entropy system is the system comprised of equations 
(2.1), (2.5), (2.6), (2.11), (2.12), (2.15), subject to the constraint u a u a = 1, 
with T a p given by (2.3), and 3? in (2.11) a given smooth invertible function. 
The unknowns to be determined are the metric g, the four-velocity u, the 
rest mass density r, and the specific entropy s. 

For previous works on the Einstein-Euler-Entropy system, the reader 
can consult the monographcs [34, 2, 8]. 

Our focus here is on the Cauchy problem, therefore we need to state 
what the initial data are. This should consist of the usual initial data for the 
Einstein equations and initial data for the matter fields. These have to satisfy 
suitable constraint equations, as we now recall. Because of our signature 
convention, the metric go on the initial three slice £ is negative-definite, a fact 
which we stress by calling the pair (S,go) a negative Riemannian manifold. 

Definition 2.4. A pre-initial data set for the Einstein equations, or pre-initial 
data set for short, is a triple {Yi,go,n), where (£,<7o) is a three-dimensional 
negative Riemannian manifold and k a symmetric two-tensor on S. A devel- 
opment of a pre-initial data set (£, go, k) is a space-time (M., g) which admits 
an isometric embedding of (E,go), with k being the second fundamental form 
of the embedding. A development (A4, g) is called an Einsteinian development 
if (A4,g) is an Einsteinian space-time, i.e., Einstein equations are satisfied 
on M. 

Let (M,g) be an Einsteinian development of (^,go,K). Then the fact 
that S is embedded into M. with second fundamental form k and Einstein 
equations are satisfied implies that the following identities hold 8 on S, 

R 90 - ML + (tr 3n k) 2 + 2 M = 0, (2.16) 

known as the Hamiltonian constraint, and 

div 90 (k - (tr 50 K)g ) + J = 0, (2.17) 

known as the momentum constraint, with /j, and J defined by 

H := JfT{n,n), (2.18) 

and 

J:=XT{n, ■), (2.19) 

where n is the unit normal of £ inside Ai, T is the stress-energy tensor, and 
T(n, •) is viewed as a one-form on S. In the above, R go , | • \ go , tr go , div go 
are, respectively, the scalar curvature, the pointwise norm, the trace, and the 
divergence, all with respect to the metric go- Notice that definition 2.4 and 
equations (2.16), (2.17), (2.18), (2.19) are general, in the sense that they do 



s See e.g. [30] for a proof. The plus sign on fi is due to our signature convention; some 
authors define J with a negative sign. 
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not assume that T is the stress-energy tensor of a perfect fluid. When T has 
the form (2.3), then (2.18) and (2.19) become 

li = .Jf(p + Q)(l-\7r(u)\ 2 J~jr P , (2.20) 



an 



d 9 



J = Jff{p + q)^J\- \ir g (u)\ 2 gQ n g (u), 



(2.21) 



where ix g : TM\^ — > TE is the orthogonal projection onto the tangent bundle 
of E. 

From (2.20) and (2.21) it is seen that, additionally to (E,go, k), to solve 
the Einstcin-Eulcr-Entropy system, one has to prescribe g, p, and a vector 
field v (which will satisfy v = ir(u) once solutions are obtained). However, as 
we adopt the point of view that the independent thermodynamic variables 
are s and r, we do not prescribe p and g directly. Rather, r and s are given 
as initial data, and then g and p are determined on E by (2.11) and (2.12), 
respectively. 

Definition 2.5. An initial data set for the Einstein-Euler-Entropy system is 
a 7-uple (E, go, k, tq, Q)j v, £&), where (E, go) is a three-dimensional negative 
Riemannian manifold endowed with a symmetric two-tensor k; r$ and <ro are 
non- negative real valued functions 10 ro,So : E — > R + ; #:MxR-> R + is a 
smooth invertible function as in (2.11); and v is a vector field on E; such that 
the constraint equations (2.16) and (2.17) arc satisfied, with fi and J given 

by 

fi = Jf(p + g ) {I - \v\ 2 go ) - Jt po, 

and 



J = ,J€(po + go)yfl - \v\ 2 go v, 



where g = g(r ,%) and p = Po(ro,So) arc given by (2.11) and (2.12), 
respectively. 

In order to state the main result, we need to impose further conditions 
that are either of physical nature or are required in order to apply certain 
analytic techniques. 

One quantity of physical relevance in non-relativistic fluid dynamics is 
the speed of acoustic waves. For relativistic fluids, we define the sound speed 
v by 

2 ( dp\ r dp . 

\dgj s p+gdr' 

where, as usual in Thermodynamics, (gj-) means treating y as a constant 
when taking derivatives with respect to x. v 2 is well defined in that, for 



9 As usual we identify vectors and co-vectors, since J is given as a vector in (2.21) but as 
a one form in (2.19). 

10 We write ?o instead of so for the entropy at time zero to avoid confusion with the quantity 
s a introduced below. 
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physically relevant equations of state, pressure cannot decrease as g increases. 
The fact that (2.22) is correctly interpreted as the speed of acoustic waves 
follows from analyzing the characteristic manifolds of the system formed by 
equations (2.5), (2.6), (2.15) and (2.11) with a given background metric; we 
refer the reader to [2, 34] for details. 

We shall require that the following inequalities hold: 

v 2 > 0, (2.23a) 

v 2 < 1, (2.23b) 

> 0. (2.23c) 

r 

(2.23a) and (2.23b) assert the physical requirements that the sound speed 
is positive and does not exceed the speed of light 11 . The interpretation of 
(2.23c) is that specific enthalpy of the fluid, defined as is positive. In 
view of (2.22), it should be noticed that (2.23a) together with (2.23c) excludes 
the possibility of pressure-free matter. Our techniques can be easily adapted 
to allow such nonetheless, leading to a theorem similar to 2.6 below 

(see section 5). 

When the initial slice E is compact and the functions involved continu- 
ous, inequalities (2.23a) and (2.23c) automatically imply that on E 

v 2 > c, and > c, (2.24) 

r 

for some constant c > 0. These conditions are needed to apply certain exis- 
tence results for first order symmetric hyperbolic systems. When E is non- 
compact, our techniques still apply, but without a bound of the form (2.24) 
and further assumptions on the initial data, the size of the time interval 
where the solution exists can tend to zero near the asymptotic region of E. 
We therefore assume (2.24), briefly commenting on more general situations 
in section 5. Also, for E non-compact, in order to accommodate conditions 
at infinity which are not too restrictive, we shall employ the uniformly lo- 
cal Sobolcv spaces 12 originally introduced by Kato [32]. Their definition 
and basic properties are recalled in the appendix A; although it comes as 
no surprise that and H s are equivalent (as Banach spaces) when E is 
compact. It is assumed that the metric employed to define the spaces _ff* ; (E) 
is equivalent, up to the sign convention, to the metric go. 

Yet another difficulty which arises in the non-compact case is that u can 
become arbitrarily close to the boundary of the light-cone in the asymptotic 
region. This would lead to a breakdown of the positivity of certain matrices 
required for our proofs, again implying that the time interval on which the 



1J It is not difficult to see that equality in (2.23b) happens if, and only if, the four-divergence 
of u vanishes, i.e., V a u a = (see e.g. [8]). When this happens, the fluid is called incom- 
pressible. Recall that incompressibility is defined in non-relativistic physics as the vanishing 
of the three-velocity, in which case sound waves travel with infinity speed. 
12 Other function spaces can be used. In particular, the "little ^ p -Sobolev spaces", £ P (H S ), 
could be employed, these being, in fact, better suited for treating more general data on 
non-compact manifolds, see e.g. [58]. 
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solutions exist cannot be made uniform. To prevent this, we impose a bound 
on the size of the initial tangent velocity v. To see that such a bound gives 
the desired control over u near the initial hypersurface, simply notice that 
in normal coordinates at a point p on the hypersurface with ^fo normal to 
E, the set of unit length future-directed vectors on T p M is given as usual 
by the upper sheet of the time-like hyperboloid which is asymptotic to the 
light-cone. 

We are now in a position to state our main result. 

Theorem 2.6. Let (E, go, k, ro, <;q, v, be an initial data set for the Einstein- 
Euler-Entropy system, with go € iJ^^E), k, v, ro, ^ E iJ.^(E), where 
s > | + 2. Suppose further that <; > 0, r > c x , > Ci, v% = 

^ 2 (To,ft)) > c\, and \ 

v \go — — c 2i for some constants c\, C2 > 0. Then there 
exists an Einsteinian development (M,g) of (E, go, n, ro, Cb; v , which is 
a perfect fluid source, with A4 diffeomorphic to [0, T E ] x E for some real 
number T E > 0. Moreover, g 6 C°([0, T E ], H^Ti)) n C 1 ([0, Te], ff^(E)) n 
C 2 ([0,T E ], J ff^" 1 (E)) ; the matter density g and the pressure p of the fluid 
source are given by (2.11) and (2.12), respectively, where the functions r and 
s m (2.11)-(2.12) belong to C°([0, T E ], n C 1 ^, T E ], H^Ql)), sat- 

isfy r > 0, s > 0, and are such that r| s = ro, s| s = ftij an d 1/2 = ^ 2 (r, s) > 0. 
Furthermore, denoting by u the unit time-like vector field of the fluid source, 
we have that u G C°([0, T E ], ^(S)) n C x ([0, Te], fl^f^E)) and tt 3 (w) = u, 
where Tr g : T^\4| s — > TE is i/ie orthogonal projection onto the tangent bundle 
o/E. 

Remark 2.7. Notice that the pointwise inequalities of theorem 2.6 all make 
sense in that s > | + 2 implies that the quantities involved are continuous. 
As discussed above, when E is compact, the hypotheses involving c\ can be 
relaxed by assuming that those quantities are simply greater than zero, 
can be replaced by the ordinary Sobolev spaces and the bound on \v\ go is 
automatically satisfied (recall that ga is negative definite). 



3. The Einstein-Euler-Entropy system in the frame formalism. 

In this section, wc shall use the so-called frame formalism to write a different 
set of equations for the Einstein-Euler-Entropy system. It will be shown in 
section 4 that solutions to the new system imply existence of solutions to 
the original Einstein-Euler-Entropy equations. As in most analytic investi- 
gations of the Einstein equations, the point of view is essentially local as a 
consequence of the phenomenon of finite propagation speed. Thus a chart U 
should be implicitly understood whenever coordinates are involved. If E is a 
space-like three surface, or a slice for short, then we shall use a slight abuse 
of notation and still denote by E the set E n U. 

In the frame formalism, the information about the metric is encoded in 
an orthonormal frame {e^j^g, where the frame is related to the coordinate 



14 



Disconzi 



basis by 

e„ = e^. (3.1) 

Notation 3.1. From now on, unless otherwise specified, all tensor fields will 
be expressed in the orthonormal frame {e M }^ =0 , with Greek letters used 
to indicate the components of such fields in the basis {e^} 3 ^. A bar over 
an index, e.g., a, indicates that it can take only the values 1,2 or 3, with 
summation of barred indices being only over 1, 2, 3 as well. The only exception 
for our choice of basis will be for the frame itself, which will be given in terms 
of the frame coefficients e A „ in (3.1). As in (3.1), capital Latin letters range 
from to 3 and will be used to denote components with respect to the 
coordinate basis; a bar, e.g., A, indicates restriction to 1,2 or 3. Indices are 
still raised and lowered with the metric g, as usual. 

By construction, in terms of the orthonormal frame {e /J }^_ , the metric 
g is always represented by the Minkovski metric g a p = diag(l, —1,-1,-1). 
The relation to the metric in the basis {gfxl^o * s gi ven by 

g AB = e\e%g^. (3.2) 

The connection coefficients T^ 7 ^ are defined via 

V Q e,g = r Q 7 (3 e 7 . (3.3) 

The condition that V is compatible with the metric then takes the form 

r/^7+ r / 7 w = a ( 3 - 4 ) 

One of the key ingredients of the formalism we shall employ is to treat the 
connection coefficients as unknowns in their own right. As T a 7 o is given in 
terms of first derivatives of the metric, treating them as independent vari- 
ables allows us to express the original Einstcin-Eulcr-Entropy system, which 
involves second derivatives of g, as a first order system 13 . 

Symmetry (3.4) will be assumed throughout. In other words, we shall 
only need evolution equations for 24 independent T's 14 with the remaining 
components explicitly defined via (3.4). 

Since it is not possible to decide from the connection coefficients alone 
whether V is the Levi-Civita connection of g [49], further conditions will be 
necessary. The requirement that the connection is torsion-free, along with the 
(once contracted) Bianchi identities and the standard decomposition of the 



13 This argument only presents the basic heuristic intuition of the method. In fact, although 
the metric will not be one of the basic unknowns, our system would be third order in the 
metric if written in terms of it. This is because the system will involve first derivatives of 
the Weyl tensor, see equations (3.9), (3.10) and (3.14e). It should be noticed that this is 
not an oddity of our formalism. If one tries to approach the Einstein-Euler-Entropy system 
in the usual formalism of second order equations, it cannot be directly solved as it stands; 
a quasi-diagonalization process has then to be carried out, leading to a system that is also 
third order in the metric [10]. 

14 Duc to (3.4), there are — independent components of F a 7 ^ in n space-time di- 

mensions, n = 4 in our case. 
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Riemann curvature tensor in terms of the Weyl and Schouten tensors, will 
be imposed as further equations of motion of the system. In this regard, we 
introduce the tensor d a ^ s (d for "decomposition" ) defined as 

d a p-ts '■= R a p~fS ~ W a p lS - g a [jS^p + g^S s ^ a , (3.5) 

where [a/3] means that the indices are anti-symmetrized; W a g s is the Weyl 
tensor; S a p is the Schouten tensor, given by 

S a p := R a p — (3.6) 
and R a g s is the Riemann tensor, which can be written as 

(O.I) 



By construction, d a s s possesses the usual symmetries of the Riemann ten- 
sor, since such symmetries are shared by W a j3lS and —g a + g^S^ a . 
Recall that the torsion of the connection is the tensor T defined via 

T a = -[«a,ep] + {Tf g - r/Je p , 

where [-, ■] is the usual commutator of two vector fields. Notice that 

T a % = -T p \. (3.8) 

Aiming at the Bianchi identities, we define the Friedrich tensor by 

F ° p 7 8 '■= W a p jS - 9 a [ T ^5]/3> ( 3 - 9 ) 

and let 

F aM := V,F\ pr (3.10) 

To understand the role of F a p 7 , consider for simplicity the case of vacuum. 
Then the once-contracted Bianchi identity 

combined with Einstein equations (2.1) yields F a p 7 = whenever the de- 
composition of the Riemann tensor holds, i.e., d a g s = 0. 
The last definition we need to introduce is 

q a ■= (p + g^V^Ua - v 2 (p + QjUaV^ - V a p. (3.11) 

q a will replace (2.6) in the new set of equations. To motivate this, suppose 
that we have a solution to the Einstein-Euler-Entropy system. Then differ- 
entiating (2.11) in the direction of it, using (2.12), (2.8), (2.15), and (2.22) 
yield 

u a W a p = -{p+ Q)v 2 V a u a . 

Using this into (2.6) then gives q a = 0. 

Tracing Einstein equations (2.1) gives 

R = -XT, (3.12) 
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where T is the trace of the stress-energy tensor. Using (2.1) and (3.12) into 
(3.6) produces 

S a fs = J(f{T a [3 - -Tg Q(3 ), (3.13) 

which is an equivalent way of writing the Einstein equations. 

We can now define the Einstein- Euler- Entropy system in the frame for- 
malism as the system comprised of 

equations (2.5), (2.8), (2.11), (2.12), (2.15), (2.22), (3.13), (3.14a) 

q a = 0, (3.14b) 

T a % = 0, (3.14c) 

d a p jS = 0, (3.14d) 

F aPl = 0, (3.14e) 

subject to the constraint u a u a = 1, with T a p given by (2.3), and where q a , 
T a ^ g, d a p lS and F al s 1 are given by (3.11), (3.8), (3.5), and (3.10), respec- 
tively. The unknowns to be determined are the frame coefficients e A a , the 
connection coefficients r ct 7 /3 , the Weyl tensor W a 3 s , the four-velocity u, 
the rest mass density r, the specific entropy s, and the matter density q. 
The usual symmetries of W a B s , and those of r Q 7 g determined by (3.4), arc 
explicitly assumed, therefore the system is written for only the 20 indepen- 
dent components of the Weyl tensor 15 and the 24 independent connection 
coefficients 16 . 

Remark 3.2. Obviously, when attempting to solve (3.14), it is not yet known 
that W is the Weyl tensor. In particular, the trace-free condition must be 
demonstrated. 

One way to motivate the choice of the Weyl tensor as one of the un- 
knowns is to consider the simpler case of vacuum. Then, the validity of 
d a p g = gives R a o s = W a p s which, upon contraction, produces the vac- 
uum Einstein equations. Since the Bianchi identities are necessary conditions 
for the solvability of Einstein equations, we also have Fp^g = ^ aW a p $ = 
V Q i? Q ( g 7(5 = in the case of vacuum. We see in this way that the usual 
vacuum Einstein equations can be recovered from the Weyl tensor 17 . 

Equation (3.14e) is sometimes referred to as the "Bianchi equation", 
and its importance in General Relativity has been long recognized. It figures 
in the Newman-Penrose's spin formalism and has been used in the study 
of massless fields [43], including gravitational radiation [48]. Friedrich has 



15 In n space-time dimensions, W a has j^n(n + l)(n + 2)(n — 3) independent compo- 
nents. As explained in remark 3.2, the trace-free condition is not assumed in the system, 
thus there is an extra freedom of ^n(n + 1) components. 

16 A11 the symmetries and which components enter in the system are described in section 
3.2, where we also write the system in a more explicit form. 

17 Othcr choices of variables are, of course, possible. Choquet-Bruhat and York derived a 
different system — also based on the Bianchi identities — where the Riemann tensor is 
one of the unknowns of the problem [11]. 
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employed it extensively, initially to obtain energy estimates in terms of the 
Bel- Robinson tensor [24, 25, 26], and later to derive several hyperbolic reduc- 
tions for the Einstein equations, both in the vacuum case [23, 22] and coupled 
to matter [21, 28], with some of these results extended in [44]. Friedrich and 
Nagy also used (3.14e) in their study of the initial- boundary value problem 
in General Relativity [27]. Using a different point of view than Friedrich, the 
Bianchi equation was used to derive first order symmetric hyperbolic reduced 
Einstein equations that treat the Riemann tensor as one of the unknowns by 
Anderson, Choquet-Bruhat and York, both in the vacuum and matter cases 18 
[1, 11]. Finally, equation (3.14e) has also been employed (in a rather elabo- 
rated fashion though) by Christodoulou and Klainerman in their impressive 
proof of the global non-linear stability of Minkovski space. 

3.1. Gauge fixing. 

The system (3.14) is overdetcrmined. In order to obtain a reduced system, 
which is determined and hyperbolic in a precise sense, a gauge choice has to 
be made. A specific decomposition of the Weyl tensor, suitable for our gauge 
choice, will also be necessary. 
Put 

7Tq/3 '■= gaf3 — U a Up, 

so 7T Q( g is the metric induced on the space orthogonal to u, with projection 
given by 

(recall that indices are always raised with g a p)- 

Remark 3.3. Despite the terminology, in general, TT a/ 3 and ir a ^ will not agree 
with the metric and projection on the t = constant space-like slices, because 
such slices are not expected to be orthogonal to u except in some special 
cases. In particular, Tr a p\ s does not, in general, agree with go- 
Letting £ Q /3 7 5 be the totally anti-symmetric tensor, with the usual con- 
vention £0123 = +1, wc define 

Expanding out, it is easy to obtain 

(3.15) 

and 

= -2ttV (3.16) 

Recall that an observer moving with a (time-like) 4- velocity V measures 
an electric and a magnetic field given by T a ^V^ and F^V^ , respectively, 



18 By the time of the completion of this manuscript, we learned of [59], where the formula- 
tion of Anderson, Choquet-Bruhat and York is used to give yet another proof of short-time 
existence for the Einstein-Euler-Entropy system. 



18 



Disconzi 



where T is the Maxwell stress tensor and T* its dual. By analogy, it is 
customary to introduce the following. 

Definition 3.4. Let W a p.ys be the Weyl tensor and W*p lS its dual, given by 
The u- electric and u-magnetic parts of W a o s are defined by 



and 



respectively 



It follows that E a p and B a p are symmetric and trace-free. With the 
help of (3.15) and (3.16), it is not difficult to verify the following. 

Lemma 3.5. The following decompositions of the Weyl tensor and its dual 
hold: 



W a p-yS = ^(j^p[-fE S ] a - UpUfoEg] a - TTafr E$]/3 + U a U[^Eg\p^ 

W* f ,„, x =2u ! „E m „.e fl „, s - _.u„ - 4ur„S fi 



(3.17) 



Decompositions (3.17) and (3.18) allow the Weyl tensor to be eliminated 
from the system (3.14) in favor of its electric and magnetic components, 
producing yet another set of equations where E a p and B a p will be unknowns; 
see section 3.2. 

Until now, {e M }^ =0 has been an arbitrary frame with respect to which 
all tensor fields have been written. In particular, there is no relation so far be- 
tween u and the frame other than the general fact that u can be decomposed 
in this base, i.e., u = u^e^. A gauge choice in our formalism will be a specific 
choice of frame — very much in the same way that a choice of gauge in the 
coordinate formalism corresponds to a determined choice of coordinates, e.g., 
wave coordinates. 

Definition 3.6. Let (U, g, u) be a fluid source. An orthonormal frame {e M }^ =0 
in IA is called a fluid source gauge if it satisfies 

eo = u, (3.19) 

with the remaining {e^}^^ being Fermi propagated along eg. By definition, 
this means 

(V eo e«,e^) = 0, (3.20) 
for a, ft = 1, 2, 3 (recall our conventions in notation 3.1). 
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Condition (3.20) means that {e a } remains orthogonal along the time- 
flow of u. Roughly speaking, this can be understood as similar to parallel 
transport, with the important difference that {e a } does not remain parallel 
but is allowed to "rotate about time-axis given by u." 

Let t be a parameter for the flow lines of u and £ be a slice. Then there 
exists a foliation of a neighborhood of £ by leaves £{ = {t = constant}, 
which are space-like and diffeomorphic to S. Let {^-i}^^ be tangent vec- 
tors on £ associated with coordinates _ . The coordinates 
can be dragged along u to give coordinates on £ t , and since u is time- like, 
u span{^}. Setting := ^ = u then gives a basis -§^a}\ =1 = 
{j£x}a=o f° r tne space-time tangent space, with coordinates {x A }\ =0 . From 
now on, it will be assumed that whenever a fluid source gauge is employed, 
the coordinates are arranged as just described, unless stated otherwise. 

The following is a simple consequence of our choice of gauge and the 
fact that g is represented by the Minkovski metric in the frame formalism. 

Lemma 3.7. In fluid source gauge, it holds that 

p A _ sA a _ m P o._n 

= g a /3 - So a 5 p, 7r a( g = —5 & p, tTq, = 8 a So a d of) . 



3.2. The reduced system of equations. 

In this section, we investigate a reduced system for equations (3.14). As 
mentioned in the introduction, its derivation, which follows that of [22, 21, 
28], is given in the appendix B. 

In light of the various symmetries involved, only equations for some 
components of the tensors involved are needed. First let us determine them. 

Taking the inner products of (3.3) with eg and using (3.4) gives 



(8 



= 0, (3.21) 



r a % = r g ^ , (3.22) 
r a % = -r a ^, ( 3 - 23 ) 

where (•)|^_, indicates that there is no sum over repeated indices, and from 
(3.4) and (3.21) wc have 

r a = r \. (3.24) 

Identities (3.21), (3.22), (3.23) and (3.24) hold for any frame. Fixing the 
gauge allows further simplifications. From now on, we shall assume that our 
frame is a fluid source gauge, unless stated otherwise. 

From lemmas 3.5, 3.7 and the symmetries of the Weyl tensor, 

E a o (= E 0a ) = 0, B a0 (= B 0a ) = 0, (3.25) 
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E &p - WoaOjgj B a p = -W 0&fll/ e Q p ^ . (3.26) 

From (3.25) and lemmas 3.5 and 3.7, it is seen that in fluid source gauge 
the Weyl tensor is determined from E & p and B & p. It is therefore enough to 
consider equations for these quantities, whereas from lemma 3.7 and identities 
(3.21), (3.22), (3.23) and (3.24), we obtain that it suffices to have evolution 
equations for connection coefficients r a ^ , r o a , T g ^ and frame coefficients 
e A - 

To write the reduced system, the introduction of yet another variable is 
needed. Write 

s a = V a s. (3.27) 

As was done for the quantities e A g , r a 7 g , E a p, B a p, we shall treat s a as 
an unknown, with the relation (3.27) to be demonstrated after solutions are 
obtained. 

Define the operator V M which acts on fields of the orthogonal comple- 
ment of eo by 

V /i j4a 1 a 2 -a« = 7r M I/ V I ,A / 3 l/ 3 2 ... A 7r ai /3l 7r a2 /32 • • • 7r 5f 01 . (3.28) 

Then 

Although, in fluid source gauge, it then follows that 

we shall explicitly write to facilitate the comparison with appendix B. 
We can now investigate the reduced system, whose equations arc 

+ (Vo ~ r o "?> A » ~ r o V^o = (3.29a) 

W- s & B + r- x ^r- s \ + r %r 5 -% - r- s \r Q %- - £ ^s a> = o (3.2%) 

,2 



ds M U"« p +g 

r dv 2 \ dp p ( p + g ( d 2 p\ _ 2 \ p 

L+ ^dr)ds^ ° S5+ l v 2 \dg 2 \ J * 00 a 



p A /p M_ p t 1 \ P A p_ ^ p_ A p M 

^VV a =0 (3.29c) 



^ f r a % - *%(r ° a ) - ^ - v 2 [r„ %-(r " a - r a * ) 

_ P -p /J i p fp £ _ p /2 \ , ,.2p M /p p_0 \ 

1 ( dg 1 dp 



p + g \ds v 2 ds 
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Jf(g + 3p)5 & g = (3.29d) 



(i 



- 3^ A r (a ° fM E- m - 2^r x \ a E- 0)fi + 5 &fj 5^8^T- x S E„ 
1 



2# B T n - v E afi - -Jfr(p + g)(T. a 



3 

A (q (q |A| A 1 v a P 



-T p \B^\ a e^^=0 (3.29f) 

$e + (p + e)r/ = o (3.29 g ) 

<9 f s = (3.29h) 
d tSa -(T \~r a \) s ^=0 (3.29i) 
9 4 r + rr/ = 0, (3.29j) 

where it is understood that in the expressions involving partial derivatives of 
g with respect to the matter variables, g is to be replaced by V (since (2.11) 
is not a part of the above system); p and v 2 are given by (2.12) and (2.22), 
respectively; and (•)r ce ui ( 8] means anti-symmetrization of the indices a and {3 

only. The unknowns to be determined are e A s , T a r o ° a , r g E^p, B a p, 
g, r, s, s a . It is explicitly assumed that E and B are symmetric, therefore 
equations (3.29e) and (3.29f) are written only for the independent compo- 
nents of these fields, say a < /3, with the remaining components defined by 
these symmetry relations. Similarly, all the symmetries (3.21), (3.22), (3.23) 
and (3.24) and the gauge condition r o ^ & = from lemma 3.7 are assumed, 
with equations (3.29b), (3.29c) and (3.29d) written only for the indepen- 
dent components and the remaining ones being defined by these symmetry 
relations. We write the unknowns collectively as a vector 

z=(e' 4 a ,r a ^,T ° a ,T a ° ^, E a p, B a p, g,r, s, s a ). (3.30) 

Remark 3.8. In order that some components of T, E and B be defined by 
symmetry, as mentioned above, it is necessary that the initial conditions also 
obey such relations. Given arbitrary E & p\ t , 1 < a < (3 < 3, one can in 
principle always define the remaining E & p\ t _ Q by E a p = E^ a . However, in the 
case of interest, E has to be determined from an initial data set, in which case 
all its components will be given on the initial slice (see proposition 4.1 below). 
In this situation, the symmetry of E at t — has to be demonstrated — and 
only then is one allowed to the write system (3.29) solely for {E a p} a< p and 
impose symmetry relations for the remaining components. Similar statements 
hold for the other fields involving symmetries. 
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4. Proof of theorem 2.6. 
4.1. Initial data. 

In order to address the solvability of the systems (3.29), we need to provide 
suitable initial conditions. These should be determined entirely by the initial 
data for the Einstein-Euler-Entropy system (as this is the set of equations 
we are ultimately interested in) and our gauge choices. We also have to show 
that the initial data for E and B, which are naturally constructed from the 
initial data set, are indeed symmetric and trace- free. Although all of this can 
be inferred from similar works treating the vacuum and conformal vacuum 
Einstein equations, as well as their coupling to the Yang-Mills equations 
[23, 24, 25, 26, 27], an explicit proof does not seem to be available in the 
literature in the case of our system with our gauge choices. It is therefore 
useful to provide it here. 

Proposition 4.1. Let 1 = (£, go, k, ro, sq, v, ,0 s ) be an initial data set for the 
Einstein-Euler-Entropy system, with Einsteinian development (A4,g) that is 
a perfect fluid source where the Einstein-Euler-Entropy system of equations 
is satisfied. Let {e a }a=o ^ e a fluid source gauge defined on a coordinate chart 
U of M, and let z be as in (3.30). Then 

can be written in terms of quantities determined entirely by I. Furthermore, 
E\ s and are symmetric and trace- free. 

Proof. It will be useful to first express space-time quantities in terms of an 
adapted frame. Let {e a }| =1 be a frame on £ orthonormal with respect to go- 
Denote by eo the future directed unit normal (with respect to g) of S. The 
frame {e a }a=o ^ s extended to M. by parallel transport in the direction of eo- 
Let {a?" 4 }^^ be coordinates on E. These coordinates are extended to M. by 
dragging them along eo, so that they are constant on the integral curves of 
eo- Denoting by x° the parameter of such curves, we obtain that {x A } A=0 is 
a coordinate system on A4. The frame and coordinate basis are related by 

~ _~ A g 

Ca ~ 6 a dx A ' 

By construction, it holds that 

e — 0' L a — u - 

Quantities expressed in terms of e a and x A will be denoted with a tilde ~, 
and the same index convention as of notation 3.1 is assumed. In particular, 
in the frame e a , the metric go is written as 

5o a/ 3 = diag(-l, -1, -1), 

and the connection coefficients are given by 
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Notice that the space-time metric is represented by the (constant) matrix 
diag(l, — 1, — 1, —1) in both frames e a and e a . This allows us to drop the ~ 
from the metric g a p , but we still write g a /3 when we want to stress that some 
expression is written in the e a basis. (The distinction has to be maintained 
though when g is written in the coordinate basis {gfr}5i=o an( i {afr}/L=o-) 
Arguing similarly to (3.21)-(3.24), 

f 7 — f — f 7 — fl 
i o p~ l & 0"^ £~ U > 

f -f P - K - 

1 a 8 — 1 a o— K ap- 

Now we proceed to relate the tilded quantities to z. Without loss of generality, 
it can be assumed that 19 

h = h on s (4 - 2) 

Write 

o d A d 

On E, e A a~§^ A can be written in terms of eg, which in turn is determined by 
gQ. As for e° a , compute 

(e a , e a ) = -1 = (^a^p ) ^ 

= (e° a) 2 + 2e A a e° affA) + e^ a e 5 a3Ai3 , 

where we have used that goo = (eo,eo) = (gfoj gftr) = 1- But since the 
projection of u onto E is u, it holds that 

d-d 

V S = (u,e s ) = ( uB q^>Z & — A ) = e A a ff(M on S > ( 4 - 4 ) 

where u' 4 = 5 and (4.2) have been used. From (4.4), we can now get g oA on 
E. Therefore, restricting (4.3) to £ allows us to solve for for e° a | s since e a is 
space-like and unit, with the sign of e° a being unambiguously defined by our 
choice of orientation. Therefore, e° a | is also determined by the initial data. 
Because e A = 5 A , we conclude that all the functions e A a are determined 
on E by the initial data (and of course our gauge choice), the same being 
true for gAB in light of (4.4) and <?oo = 1. 

The change of basis from {e a } to {e a } is given by a Lorentz transfor- 
mation A: 

e a =A" a e /1 . (4.5) 
Taking the inner product of (4.5) with J^a and restricting to S produces 

e B a g B A = A p ofi B p9BA on S, (4.6) 



19 Then on £ it also holds that g^B = 9AB> but ^ w '" be convenient to write the ~ in the 
metric as a way to keep track of the changes of basis. 
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where (So, J=i ) = has been used. Prom our previous relations, it follows that 
all quantities on (4.6), except possibly the A's themselves, are determined by 
the initial data on E. Viewing (4.6) as a system for the hP a: (a = 0, . . . , 3, 
p, = 1, 2, 3) (which will be invertible since the matrix of the system is that of 
the change from -^~a to up to a lowering of the indices) shows that Q | s 
is entirely determined by the initial data and the gauge choices. Recalling the 
identity 

A%A^= 57 , (4.7) 

then gives 

(A° ) 2 = 1 + £(A%) 2 , 
0=1 

so that A is also determined only by the initial data when restricted to 
E. Notice that the sign of A is positive because A belongs to the proper 
Lorentz group 20 . From (4.7), the remaining A's are given in terms of hP a 
and we conclude that A| s is determined by the initial data. 

Next we investigate the connection coefficients. Taking the inner prod- 
ucts of (3.3) with eg, 

ds^a^p = (Vae^ea). (4.8) 

Writing the frames on the right hand side of the above expression in terms 
of S Q via (4.5) leads to 

T Q % = 7 ff 7 *A" a e„(A" p ) + g TtT k° 7 ^ A" ^ „f „ \. (4.9) 

On E, the coefficients T v are determined by go. Then, by virtue of (4.1), 
the construction of So and e^, and our previous relations involving A, it 
follows that all terms on the right hand side of (4.9) are determined by the 
initial data, except possibly those involving the derivatives of the Lorentz 
transformation in the direction of So. To see how such terms are determined, 
recall that in fluid source gauge T Q * ^ = 0, so that (4.9) gives 

= A 7 <7^A° e (A° ) + g B<r A° ^A? S p (A* - p ) 

Because A > 1, we see that 

A° 7 ^e (A%)| s (4.10) 

can be written in terms of quantities determined by X. Using then (4.10) into 
(4.9) with a M> a, fi n- /3 and £ H- £ shows that r_ ? ^ 

solely in terms of quantities coming from X. The coefficients r Q ^ 

r_ ^ are similarly determined: from r o ° = and (4.9) one obtains that 



is also written 
and 

s 



-"Since u belongs to the inside of the future light-cone. 
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eo(A° )| is determined by X; using this fact into the expressions for r o ^ 

and r g ^ shows that the same is true for these quantities restricted to E, 
and hence the claim follows upon evoking (3.22) and (3.24). 

The restrictions of g, s and r to £ have the desired properties in that 
they are just scalar functions. For s a , our gauge choice and (2.15) imply 
So = 0. Then, recalling that -A^ = eo, 



A ds_ 



6 ™ ~ A j 



which is determined by X on E by the above results for e ^ and noticing 
From definition 3.4, 

E af) = Woaopu ^ + W 0aS pU°W + WpcO^U + Wpcurp&iF. (4.11) 

On T, we have vf 1 = , while u° | s is computed from the w and the normaliza- 
tion condition u a u a = 1. Therefore, to show that E a g is written in terms 

of quantities determined by X, we only need to investigate the components 
of the Weyl tensor in (4.11), and by its symmetries, it suffices to do so for 
components of the form W p g, W oaff j and W ms p. 

From the Gauss equation, the decomposition of the Riemann tensor 
(i.e., (3.5) with d a ^s = 0), Einstein equations, and (3.12), we obtain 

- (3) %« - K- x X Kp- s + k x gK B - x on E, 

and 

+ - ^Tg^gtf + ^Jf(Tp p - ^Tg^g^ (4.13) 

- i^(T ffp - on S, 

where ^ Rp$ and ^R^asp are respectively the Ricci and Riemann curvature 
of (S, go)- hi a similar fashion but using now the Codazzi equation: 

WW = (3) V s ^ fi - (3) V,gK ffa + ijrf 0( gg aff - l -WT m g^ on E, (4.14) 

where is Levi-Civita connection of go. 

From (2.3), (2.11), (2.12), (4.11), (4.12), (4.13), and (4.14), we obtain 

that E a p is a tensor solely determined by X, which is symmetric and trace- 
free by the constraint equations (2.16) and (2.17). The same holds for E a p\^ 
by the invariance of the trace and the properties of v previously shown. 
By an analogous argument, a similar statement holds for B a p\y,. □ 
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Definition 4.2. By proposition 4.1, given an initial data set, a choice of (fluid 
source) gauge uniquely determines initial conditions for the reduced system. 
These initial conditions for the reduced system are henceforth called a reduced 
initial data set. 

For practical applications, e.g., to numerically solve the equations, propo- 
sition 4.1 tells us how to arrange the initial data. Given a negative three di- 
mensional Riemannian manifold (£, go), choose coordinates {-§^a and an 
orthonormal frame {ea}| =1 . Declare a metric on [0, T] x £ by g = dt 2 + go, 
where the tilde emphasizes that we identify [0, T] with the time coordinate 
a; (so that eo = gfo), and the metric is written in terms of the resulting 
coordinates on [0, T] x S. Using this metric and the (given) three- velocity v 
we determine the four velocity u on £ by the condition (u, u) = 1 (notice that 
flowing along u produces the parameter t, i.e., xo). Besides, we complete u 
to an orthonormal frame {e Q } Q=0 a ^ the instant x° = 0. The relations given 
in proposition 4.1 can now be used to define the remaining quantities such 
as A q , T a ^ p, T a ^, etc., on the initial Cauchy surface. 

4.2. Well-posedness of the reduced system. 

Short-time existence for the reduced system is a direct consequence of the 
way it has been set up, at least under our hypotheses. In fact, the gauge 
choice and construction of (3.29), originally devised by Fricdrich [21], are 
motivated exactly by the attempt of obtaining a reduced system that is sym- 
metric hyperbolic, in which case well-known results can be applied. There 
are, however, one subtlety and one observation, that have to be dealt with. 
First, the initial conditions for the e A s , T, E and B involve a different num- 
ber of derivatives of the initial metric go, hence they belong to £f^(E) with 
different values of t. The usual techniques of symmetric hyperbolic systems, 
however, yield solutions with the same regularity for all the unknowns (see 
e.g. [39, 56]), which in this case would be that of the less regular initial data, 
namely, E and B. This docs not give the desired differentiability for the frame 
coefficients, nor for the metric. Second, although at first glance the matrix 
coefficient of appears to be a diagonal matrix with entries either 1 or i/ 2 , 
the "spatial" derivatives involve derivatives in the direction of x° hence 
contributing to the zeroth matrix coefficient. 

Proposition 4.3. Let X be an initial data set for the Einstein- Euler- Entropy 
system satisfying the hypotheses of theorem 2.6. Fix a positive real number 
Te and consider the system (3.29) defined on [0,Te] x E, where x° is iden- 
tified with the time coordinate on [0, Te] ■ Let To be the reduced initial data 
set determined on {0} x £ by I. Then there exists a unique solution z to 
(3.29) on some time interval [Q,T' E ], < T' E < Te, and satisfying z\ t=0 = 
2b. Furthermore, e A a G C°([0, T' E ], H^{H)) n ^([0, T E ], #»,(!])), e° s £ 

c fl ([0 1 ^] I ffi(E))nO 1 ([0 I f B ] I ^- 1 (E)) ) r 5 %, r \, r a % e c°([o,r E ], 
^(s))nc 1 ([o,T^],F^- 1 (S)) ) E a p, b &p e c ([o,r^],ff i 7 1 (i]))nc 1 ([o,^], 
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C o ({0,T' E },H. 



s a 6 C°([0,T E ], 

,(s))nc 1 ([0,ri 



ff: I - 1 (E))nc 1 ([o ) n] ) ^- 2 (S)) ) e , r ; 

B J.^i _1 (S)). 



Proof. First, we claim that (3.29) is symmetric hyperbolic with respect to 
t on the initial hypersurface and remains so as long as v 2 > 0, and the 
slices St = {t = constant} are space-like with respect to the quadratic 
form g t induced by the frame coefficients. Notice that symmetry here means 
symmetry of the matrix coefficients M A of the derivatives • Therefore we 
have to first change basis via (3.1). Other than the first term in each equation, 
the derivative also figures in the terms involving ep in equations (3.29c), 
(3.29d), (3.29e) and (3.29f), where in these last two equations the contribution 
of ep comes from the covariant derivatives. Expressing all derivatives in (3.29) 
in terms of gives that the term in <9 t can be written symbolically as 
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(4.15) 



where fnt^ is the matrix part corresponding to 



2 C M C £ 



d t B, 



From (4.15) and (4.16), it is seen that M° = M l is symmetric. Symmetry of 
the remaining M A , A — 1,2,3, is similarly verified. 
The quadratic form g t is given by 

3 



2 C fJ. 



d t E ri 



(4.16) 



where the coefficients f a A are defined via 

From these constructions, we obtain that the characteristics of the system 
are non-zero multiples of 

where Ky, . . . , K4 are positive integers. It follows that the system is symmetric 
hyperbolic as long as v 2 > and gt remains negative definite. 
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Consider now the problem on a local patch [0,Te] x U. From our hy- 
potheses and the constructions of proposition 4.1, it follows that the initial 
data To = z(0, •) is such that 

e° a (0, •) £ H'(U), e A & (0 1 •) G H S+1 (U), (4.18) 
( Q> r,s)(0,-)eH s (U), (4.19) 

(T a %,T ° a ,T a %)(0,-)eH s (U) 7 (4.20) 

(^,%« a )(0,-)Gff'" 1 (lO- (4-21) 

e° a (0, •) is only in H s because it depends on v (see (4.3) and (4.4)). From 
(4.18)-(4.21), wc conclude that Zo| s G ff s— 1 . This is enough to apply the 
theory of quasi- linear symmetric hyperbolic systems as in [20, 32] (recall 
that s > | + 2), whose hypotheses are satisfied due to the above positive 
definiteness of M*. Shrinking U if necessary, we obtain a unique solution zjj 
in C^T^iP" 1 ^)) n C x {%T' E \,H a -\V)) for some < T' E < T E . 

In order to obtain the desired regularity, we shall use a bootstrap ar- 
gument. Consider the system for the frame coefficients formed only by equa- 
tions (3.29a), where the T's now enter as inhomogeneous or lower order terms 
given by z v , thus they are in C°{[Q,T E ], H 3 ' 1 ^)) n ^([0, T' E ], H S ~ 2 {U)). 
This is just a first order symmetric linear system for the frame coefficients, 
but there is a mismatch between the initial data, which is in H s by (4.18), 
and the lower order/inhomogeneous terms, which arc in H s . The results 
of Fischer-Marsdcn [20] deal precisely with this situation, and we obtain 
therefore a unique C°([0, T E ], H S (U)) n C^flO, Tg], H 8 ' 1 ([/)) solution 21 . By 
uniqueness, this solution agrees with that of zjj for [0,min{T^,Tg}], and 
shrinking the intervals if necessary, we can assume T' E = T E . 

Next, consider the system of equations (3.29g), (3.29h) and (3.29j). As 
before, the initial data, given by (4.19) is in H s , whereas the coefficients 
are only in Notice that although the matrix coefficient of is the 

identity, depending on the equation of state, this system is semi-linear due 
to the presence of the pressure in (3.29g) (see (2.12) and (2.11)). In any case, 
the results of [20] still apply, and we obtain that g, s and r are also in H s . 

A similar argument can be applied to the system (3.29') comprised 
of equations (3.29b), (3.29c) and (3.29d). The iJ 3 " 1 terms E & p, B a p, s„, 
and the now H s terms g, s and r given by zjj enter in the system (3.29') 
as inhomogeneous or lower order terms, whereas the initial data for (3.29') 
given by Iq is in H s by (4.20). We write the system once more in terms 
of the derivatives g§x, obtaining a semi- linear system where the coefficient 
matrix N* of involves the frame coefficients e A & and the sound speed 
v 2 . which is given in terms of s and r by (2.22). From the above argu- 
ments, we obtain that N* is positive definite and is in H s . Evoking the 
results of [20] one more time gives that r a ^ T Q a , r a °^ are in fact in 
C°{[0,T' E ],H s (U)) nC\[0,T E },H s - 1 (U)). Using this improved regularity of 



21 The results of [20] apply to a large class of quasi-linear equations. Allowing for less 
regular lower order terms was, as the authors acknowledge, one of the goals of the paper. 
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the connection coefficients again with (3.29a) and (4.18), for A = 1,2,3, 
finally gives e A a € C°([0,Tg], H S+1 (U)) n ^([OjTg], H S (U)). 

We now obtain the result on [0,T^] x E by a standard gluing proce- 
dure. Uniqueness guarantees that solutions constructed from different patches 
£/and U' agree on the domain of dependence 22 of U D U'. The time interval 
[0, T' E ] can be made uniform due to the uniform conditions on the initial data, 
and the local in time, global in space, solution z will belong to the desired 
H*, spaces by the way these spaces are constructed out of the Sobolev spaces 
of maps defined on local patches. □ 

The geometric meaning behind the definiteness of the matrix M is easy 
to grasp. If u were hypersurface orthogonal, then e a would be tangent to E t , 
and M* would be a diagonal matrix with positive entries. By continuity, we 
would expect M to remain positive definite as long as u is sufficiently inside 
the light-cone. 

We also notice that the above bootstrap argument for the regularity of 
some of the components of z works because of the particular form of system 
(3.29), which can broken in several sub-systems that are "mildly coupled" 
among themselves. Since several of the quantities involved have direct phys- 
ical meaning, it would be interesting to see if such split into sub-systems 
can have an useful physical interpretation, perhaps in terms of some effective 
notion of weak coupling among certain quantities. 

Corollary 4.4. The solutions E and B constructed in proposition 4-3 are 
trace- free, and g satisfies (2.11). 

Proof. Tracing equations (3.29e) and (3.29f), we obtain a first order symmet- 
ric hyperbolic system for the traces of E and B. Since -E Q a | s = = B a a | s 
by proposition 4.1 and (3.25), by uniqueness these tensors remain traceless. 

Locally s and r are written in terms of the coordinates s = s(x , ... , a; 3 ), 
r = r(x°, ... ,x 3 ). From (3.29g), (3.29h), (3.29j), (2.12), we obtain 

-&(r,s) = d tg , 

which implies g — £P(r, s) since this holds at t — 0. □ 
4.3. Propagation of the gauge. 

Letting e A = S A and A4 = [0,Tg] x S, from proposition 4.3, we obtain a 
space-time (A4,g) with the metric given by (3.2). Notice that g agrees with 
go on E because of (4.17). (M,g) is turned into a fluid source by setting 
u = eQ. The components of the fields that are not given in proposition 4.3, 
e.g. E a Q etc, are defined by their corresponding expressions in fluid source 
gauge, e.g. (3.25), and their symmetry relations. The frame {e Q }^,_ is then 
a fluid source gauge, with the quantities of proposition 4.3 being exactly field 
components written in this gauge. All other quantities throughout this section 
will be written with respect to this frame unless stated differently. Moreover, 
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we shall also assume the hypotheses of theorem 2.6, so that the results of the 
previous section will also be used throughout. The coordinates are arranged as 
explained after definition 3.6; in particular eg = Jj. By construction, (3.4) is 
satisfied and the connection associated with g is compatible with the metric, 
but it is not known at this point whether it is torsion free. In particular, in 
all expressions below involving a covariant derivative, it is to be understood 
that V is such a connection and not the Levi-Civita one, at least until the 
torsion free condition is demonstrated. From E and B, we define W and W* 
by (3.17) and (3.18). W has then the usual symmetries of Weyl tensor and is 
trace-free by corollary 4.4, but it is not yet known that W is the Weyl tensor 
of the metric g. With u and g known and p given by (2.12), we define T a p 
by (2.3). Recall (3.7) and define d a ^ s via (3.5), where S a p in (3.5) is given 
by (3.13). Define also T a 7 /3 , F a 0jS , F a01 and q a , by (3.8), (3.9), (3.10) and 
(3.11), respectively, where v in q a is given by (2.22), with s, r, and g being 
those of proposition 4.3, which satisfy g = S?(r, s) by corollary 4.4. We now 
proceed to show that (3.14b), (3.14c), (3.14d) and (3.14e) arc satisfied. In 
order to do so, we shall derive a symmetric hyperbolic system of equations 
for these quantities and show that they vanish on the initial slice. 

Remark 4.5. At the risk of being repetitive, we stress again that when refer- 
ring to equations such as (3.17) and (3.13), it should be understood that they 
are being used to formally define W, S a p etc, from the quantities obtained 
from proposition 4.1. Notice also that because we do not know that V is 
the Levi-Citiva connection, the torsion tensor will have to appear in several 
manipulations below. 

Lemma 4.6. With the above definitions, d a p s enjoys all the symmetries of 
the Weyl tensor, T a 7 a = —T f3 1 a , and (2.5), (2.15) and (2.8) are satisfied. 
Furthermore, qq = 0, and 

0, (4.22) 
0, (4.23) 
q a . (4.24) 



T p - 

'0 a ~ 

d & - = 

U /307 

V^T — 



Proof. The symmetries of d a ^ s and anti-symmetry of the torsion tensor are 
direct consequences of their definitions and the fact that W a j3lS has these 
symmetries. In fluid source gauge, (4.22) and (4.23) are equivalent to (3.29a) 
and (3.29b), respectively, while (2.5), (2.15) and (2.8) are the same as (3.29g), 
(3.29h) and (3.29j), respectively. 

Since the pressure is a function of r and s by (2.11) and (2.12), differ- 
entiating p with respect to t and using (2.22), (3.29h) and (3.29j) yields 

d t p+{p + B yV/ = Q. (4.25) 

(4.25) is equivalent to qg — in our gauge. Computing V^T^q, from (2.3) and 
using (3.29g) (or cquivalcntly (2.5)), 

V^T MQ = (p + e)«"V^ Q + u^V^p - V a p, (4.26) 
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which in light of (4.25) and our gauge conditions, produces (4.24). □ 

The next lemma and the proposition that follows will be the main in- 
gredients in proving the propagation of the gauge. Although both proofs are 
heavily computational, they follow the same lines of [21, 28, 27]. 

Lemma 4.7. The following relations hold: 

d t F Qa0 - \e & ^V,F 0X , + df [v " a W a]li&0 - W^d*^ 
~ 2% ^Vo-W^gg + d? [o^ajfi^w + d^ J ^S ^ + T^ F oa o 
- 7 r /o^oAo(Voe P a o - \^e% & )e/ s F^ s = (4.27a) 



+ 2^ 

1 niR „ 1 



[a fl/^ ^ Oo-A ^ 

1 

4"" ^ 1 ^" °7 " ao/3 



^(p + ^ s V ar (4.27b) 



(0/3 7 ) 

- E r o Va % " E "o EVo^A = (4.27c) 

<^7> <70> (0£> 

9 *^7 + E Rfi f>ooTp % + \e m \F- XTi + 9- XT q,)e^ + T o ^d% - 

(0/37) 

- E V^o ~ E r 7 = ( 4 - 27d ) 

<P70> (P0/3) 

(0| 7 ) 

+ g Xrqi )e PT « + <^r V* + ^ % v P CT?0 + s,% v p 

+ ^5^r/Xoo^=0 (4.27e) 
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Se^ [ ]T R° , aa T- % + \s &m \F- XTi + 9xT q,)e v 

(a/37) 

lr d a 



0t£ 



6 



^cr 1 XI 7CT 

1^7 ^ 1 p a a PT-a " 1 " 7 cr vafl 



Vr'd - V x 

/ j ol v (j^7 / , ^ /3 u <J~ta 



(a/37) 

1 1 a£ 



(f7Q) 



rs-tf^* - V r * d° -3 

/ j 7 u caff 

(pa/3) 



9 r <T A o 7 a ^a - -7r a M 7r / g ,y (d A opp + d A pp0 + d x mp )q\ 

+ 2^ 2 (p + e )r a d p 0igp - 2^ 2 (p + e)r o y oap - (^ + a te )d a0 ^ 
+ v\ v + g )(T p \ - r a ^K oX _ + (p + e )r a A d° - (p + g )r p \J 

+ e a (V (p + g)) - (V (p + g)) d** mfl 



Aoa 







(a ^%)A, " «(a(^ , _) + iw^dO ^(^) 



(4.27f) 



1 1 



2p+g 



" Q p, 0A£ x ffOA 



i /3 A" OmC P A" flOa) + 2 e ( a 2^ « P) CT 0',3 



'/3 7 



(O/37) 



+ 



y { J"e^s m \F Xr(+9XTq ,) + \e { ^(d% yM T \ 



- A ^- < i - 2/ ~ 1 -0, 



+ (r a A o - r A a )^ - 2^ 2 (p + )cT 



Oa-y 



(4.27g) 
(4.27h) 



where (-)(q|^|^) (Vesp. (-)[a|^|^] ) means symmetrization (resp. anti-symmetrization) 
of the indices a and f3 only, Yltafii) indicates sum over cyclic permutations 
ofaPj, and 



(a/3) 

T T x a A X p 1 := T T x a A X p 1 + T T x l3 A aXl + T T x 7 A a p X . 

(a/3-y) 
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Proof. We start computing V^-F^q,^. Commuting the covariant derivatives 
and using the symmetries of W 11 a ^ and S a p we find 

— ^ p a^pv+R a pS^u - R 11 u^aSpLi + R v 0^afi 

Using (3.5), (3.13) and (4.24) and the various symmetries of the tensors 
involved, the above becomes (recall that T is the trace of T a p) 

w, a , 3 = - dV'^w + W ^JT + \% w 



(4.28) 

On the other hand, using (4.24) into F^ a p gives, after some contractions, 

-\jTq a = ^V^F^y^ - tt^F^uW. (4.29) 

Then, from (3.17), (3.29e), (3.29f), (3.29g) and (4.29), we obtain, after some 
algebra, 

F^p = r K a X F a \ u u (T u v up - ■K li x F a \ v u y u v u a + --K tia ^p X F tl \ l ,u^u u 

- \-K^ a X F^ v u»u» - ^u^ a0 e^n a \ T ^F^ (4.30) 



CTT 



Computing V M F Ma ^ from (4.30), using the resulting expression into (4.28), 
recalling our gauge conditions, and evoking lemma 4.6 leads to (4.27a) and 
(4.27b) after suitably choosing the indices to correspond to the ones of those 
expressions. 

When the torsion of V does not necessarily vanish, the first Bianchi 
identity takes the form 

E = E ( v «^ % r a v 7 \) ■ (4-3i) 

(a/37) (a/37) 

After evoking (3.17), using symmetries, and setting a = 0, /3 = /J and 7 = 7, 
this expression simplifies to 

VoT^ % = -v t, " - v,r + £ (>o^ + r A ,-r, <\ 

(007) 

In view of (4.22), this gives (4.27c). 



34 



Disconzi 



Taking torsion into account again, the second Bianchi identity reads 

E V Q i?% 7 + ]T R\ Xa % \ = 0. (4.32) 
(a/37) (a/37) 
From the symmetries of W a ^ lS , d a p 1& and (4.31), we have the identity 

E V o n^ + ^»WV T r fA ,=0. (4.33) 

(a/97) 

From (4.32), (3.5), (3.17), (3.18), (4.33), (3.13), and (4.24), we get, after some 
simplifications, 

(a/97) (a/37) (4.34) 

+ ^ Q/37 A (^+.9Ar^)e^ re = 0. 

Setting a = 0, fj, = p, v = D, /3 = /3 and 7 = 7 in (4.34) and using (4.23) 
produces (4.27d). 

From the definition of q a , compute V [5(7^1, use (3.29d) and v 2 > to 

find 

(p + e)rf a o^- V [a9/3]- (4-35) 

which implies 

Then (4.36), (3.15) and the symmetries of d a p s imply the identity 

1 
2 

From (4.34) with a = 0, fj, = u = v, /3 = j3 and 7 = 7, 

eo(d° QM ) + e ? (d° p , ) + e,(d° mp ) + £ R\ M T f3 % + V Q \d% p . 

(0/37) 

+ ^p"a d ' J p70+ T j ad Cr ^ p~ E r CT 5 rf0 <T^7 _ E ^ <n"0 

( P( 9 7 > <P70> 

- E F 7 + ^ HFxrt + g-xr*)^ = 0- (4-38) 

pop) 



9^^(d° m p)-9 W e (dP m ) = s ^e x (e^d° (4.37) 



A' 



Using (3.21), (3.22), (3.23) and (3.24), observing that R° P(tA = R p 0i 

d° i>a\ = ^Oo-A' usm g (4.37) into (4.38), and contracting in v and 7 pro- 
duces (4.27e). 

From (3.29c) and the definition of q a we find 

(p + g yd\ aiy = \7 [Q q a] . (4.39) 

From (4.35) and (4.39) we then obtain 

eo(V [5 ^]) = -d° &0 pe Q (p + g) - {p + g)e (d° &o0 ), (4.40) 

V Q V [0 ^] =e 5 (V [0 ^ ] ) -r a x V [o q x] -T & x Q V [m , (4.41) 
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V^V [ofc] = ep(V [0 q &] ) - \V [0 q x] - A V [Afe] . (4.42) 

Recall the following identity for the Lie derivative C 

Ceo^laQp] = VoV [a ^j + V^rg M Q + V^^r- (4.43) 

Now compute the left hand side of (4.43) directly from (3.11), use (4.40), 
(4.41), (4.42), contract the resulting expression with e p a ^ , and evoke (4.34) 
once again to find (4.27f). 

Next, set a = 0, fi = 0, v = v, /3 = j3 and 7 = 7 in (4.34), contract with 

e a 1 symmetryze on a and i>, and use (4.35) and (4.39) to obtain (4.27g). 
Finally, (4.27h) follows from (4.39) and lemma 4.6. □ 

Next, we show that the tensors d a g s , F a ^ and T a 1 p vanish on Ai. 

Proposition 4.8. With the above definitions, 

d a Sj6 = 0, F aPl = 0, T a 1 B = 0, and q a = 

on M. 

Proof. In light of our gauge choice, lemma 4.6, and the symmetries involved, 
several components of the above tensors vanish identically on Ai. Taking into 
account the symmetries of the remaining components, it is seen that to show 
the proposition it suffices to prove that the components 

T & 7 ' 0, Foao, F m p, F a ^, d° a0( g, d° a p-, d^-^., and q & (4.44) 

vanish on M. 

Recalling that F a p 7 = —F ai p, equations (4.27) can be viewed as a 
system for the quantities 

T- 7 p, F oa o, F m p, d? .p., d v mD , £p Q/3 rf° a0 ^, e (v^d° a)M' and q& - ( 4 -45) 

Arguing as in the proof of proposition 4.3, we obtain that (4.27) is a first order 
symmetric hyperbolic system. Uniqueness of solutions then implies that all 
quantities in (4.45) vanish on M if they vanish on S. To show that this is the 
case, it is useful to introduce the adapted coordinates {x A }\=o and lrame 
{e/i}^=o as m tnc P r001 01 proposition 4.1, and we shall employ the same 
notation and conventions as used there 23 . 
From (4.22), it follows that 



Using (4.5), this gives 



r a >« = o. 



% v°o + r a v a o = °- 



But since the connection on three manifold (£,go) is torsion free, we have 



T 1 - 

'a 



= 0, whereas 71 s 

5 a P 



by (4.1). Therefore, T Q 



since 



A° 7^ 0. We conclude that T a 7 a, and thus T a 7 g , vanishes on £ 



23 As some identities of the proof of proposition 4.1 will be evoked as well, it should be 
noticed that these are still valid assuming only the hypotheses of this section. 
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From (4.30) and the way the initial data was constructed in proposition 
(4.1), we find, after a somewhat lengthy but not difficult calculation, that 

F oao = = F mfj on S. (4.46) 

Using corollary 4.4 and lemma 3.5, we obtain, with the help of (4.24), 

In light of (4.46) and (4.30), this gives 

q~ t = on S, 

which then implies, upon employing (4.30) one more time, that 

F &M = on S, 

where (4.46) has been used. 

Finally, using the Gauss equation and arguing as in proposition 4.1, we 
conclude that the constraints implied by the decomposition of the Riemann 
tensor (i.e., (3.5) with d a ^ $ = 0) arc satisfied on S. This, combined with 
(4.23), gives the vanishing of d a g § on E. 

We conclude that the quantities (4.45) vanish on Ai. This implies, evok- 
ing (4.30) once more, that F & ^ = also holds on Ai. The remaining com- 
ponents in (4.44) also vanish due to the identities 

which are verified by inspection with the help of (4.36), (3.15) and (3.16). □ 

Remark 4.9. Notice that, as in the familiar case of wave coordinates, propa- 
gation of the gauge requires that the constraint equations be satisfied. 

We already know that the connection associated with T^ 7 ^ is metric. 
Since it is also torsion-free by proposition 4.8, we obtain: 

Corollary 4.10. The connection defined by the coefficients T^ 7 ^ is the Levi- 
Civita connection of the metric g. 

4.4. Solution to the original system. 

Showing that the solution z of the reduced system yields a solution to the 
original Einstcin-Eulcr-Entropy system is now a matter of unwrapping all 
our definitions. 

Proposition 4.11. Let W a g~ s and S a p be as in section 4-3 and g the metric 
constructed out of the solution of the reduced system given in proposition 4-3. 
Then W a s ^ s and S a p are, respectively, the Weyl and the Schouten tensor of 
the metric g. Furthermore, equations (3.14) are satisfied, and the quantity s a 
from proposition 4-3 is in fact the derivative of s. 
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Proof. The connection given by T a 7 „ is the Levi-Civita connection of g by 

corollary 4.10. Thus, if we denote by W a p s and S a p the Weyl and the 
Schouten tensor of the metric g, we have 

R a p-yS = W a p lS + g a [yS^p - g^Y, S5] " 

which in turn equals W a ^ s + g a [y S s] ^ - g^S s] a since d a ^ s = by propo- 
sition 4.8. Hence, tracing the equality 

W a M5 + g a b S s]/3 - g f3h S s] a = W a M5 + g a h S 5]p - g^S s] a (4.47) 

and using that W is traceless in light of corollary 4.4, we obtain that S a p is 
indeed the Schouten tensor, which then implies W a p s = W a ^ s by using 
(4.47) again. Then (3.13), with S a p being the Schouten tensor, also holds. 
The remaining equations of (3.14) are satisfied by propositions 4.2 and 4.8, 
and corollary 4.4. Notice that these results show the validity of the (3.14) in 
fluid source gauge, but by the tensorial nature of the equations, they hold in 
any frame. 

Put Sq = V a s. Then, since u a V ' a s = by (3.29h), we obtain £„s ce = 0, 
where C is the Lie derivative, which in turn implies 

d t s a -(T \-T a %)s^ = 0. 

From (3.29i) and the construction of the initial data (proposition 4.1), it 
follows that Sq, = s a . □ 

Proof of theorem 2.6: By proposition 4.11 and (4.24) we obtain that the 
equations of the Einstein-Euler-Entropy system are satisfied. Notice that 
u a u a = 1 holds by the way the solution was constructed; so V M r MQ = 
in fact implies (2.5) and (2.6). That M. ~ [0, Tg] x S is indeed an Einsteinian 
development follows from the fact that the constraint equations are satisfied, 
and by construction (M.,g,u) is a perfect fluid source. 

Let {e^^-g be a perfect fluid source gauge, with the coordinates {^ A }^ = o 
arranged as explained below definition 3.6. Since e A a is in C°([0, Te], H^ 1 ^)) 
nC* 1 ^,!^], #*,(£)) and e A Q = 5 A , by (3.2) wc conclude that the inverse 
of the one-parameter family of metrics gt induced on E f = {t = constant} 
belongs to C°{[0,T E ],HZl l (E)) n C x ([0, T E ], H^(E)), and so does g t itself 
because s > § + 2. Since d a M5 = and E & p, B &j3 e C°([0, T E ), H^ 1 ^)) n 
C 1 ([0,T E ],H^(S))nC a {[0,T E ],H^ 1 (S)), by (3.17), (3.25), (3.26), propo- 
sitions 4.8 and 4.3, we conclude that g G C°([0, T E ], lf^ +1 (S)) n C l ([Q,T E ], 
H^(E)) n C^Te],^ 1 ^)). The four-velocity u belongs to C°([0,T E ], 
iJ^(E)) n C 1 ([0, Te], jy*r a (E)) and has the correct projection because of 
(4.4) and proposition 4.3. 

The functions r and s have the desired regularity and take the correct 
initial values by the way they have been constructed, and g and p arc given by 
(2.11) and (2.12) as a consequence of corollary 4.4, which implies that v 2 also 
has the correct form. By continuity on the time variable and the hypotheses 
of the theorem, we obtain that r > and v 2 > for small T E . We cannot 
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have s(j>) < for a point p near the initial Cauchy surface because this would 
contradict u a V a s = and <?o > 0, hence s > 0. □. 

Remark 4.12. The less regular frame coefficients e° g do not affect the reg- 
ularity of the space-time metric because they contribute only to the mixed 
entries g ^, which are gauge terms therefore having no direct physical or 
geometrical meaning. 

5. Further remarks. 

The case of barotropic fluids is treated as a particular case of theorem 2.6, 
at least as long as (2.23a) and (2.23c) hold. Although the condition v 2 > is 
violated by pressure- free matter, in this situation, equations (3.14) simplify 
considerably; a reduced system which does not require the introduction of v 2 
can be derived [21] , and a system for the propagation of the gauge, which does 
not involve v 2 either can also be constructed [28]. The arguments of section 
4 can then be reproduced, yielding a statement analogous to theorem 2.6 for 
pressure-free matter. Notice also that when a fluid is iscntropic, the entropy 
can be treated as a parameter in the equation of state, and the equations of 
motion take the form of those of a barotropic fluid. 

It should also be noticed that condition (2.23b) has never been used. In 
fact, such inequality is necessary due to causality, but it plays no role on the 
well-posedness of the Einstein-Euler-Entropy system. There is at least one 
situation where it may be desirable to consider equations of state where v 2 < 
1 is not satisfied, namely, the construction of appropriate gauge conditions 
for the vacuum Einstein equations. In this situation, setting Jf = in our 
system, the Euler equations decouple, and the role of the four-velocity u is 
to fix the gauge. As Fricdrich has pointed out [21], this procedure may be 
particularly important in numerical treatments of Einstein equations, where 
ever more sophisticated gauge choices are crucial for accurate results. 

Finally, wc remark that although theorem 2.6 does not provide an ex- 
istence result for the fluid body discussed in the introduction, when one is 
interested solely in its behavior near a small compact set O C S, the hy- 
pothesis that ro is uniformly bounded away from zero can be replaced by the 
condition that ro is positive in D, and decays sufficiently fast on its comple- 
ment. As mentioned in section 2, this will not generally yield a uniform time 
span for the solution, but a uniform Tg > will exist in the neighborhood 
of n. Whether this suffices for studying the dynamics of the fluid near il will 
obviously depend on the particular application one has in mind. 

Appendix A. The uniformly local Sobolev spaces. 

Below we review the uniformly local Sobolev spaces originally introduced by 
Kato [32] . Some of their properties can also be found in [8] . 

Let (M, 7) be a Riemannian manifold. For any open set U C M, let 
H S (U, 7) be the Sobolev space of tensors of a given rank defined on U, with 
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the derivative V and the measured of integration [i used to define the Sobolev 
norm being those of the metric 7. Recall that -ff ; s oc (M, 7) is defined as the 
space of tensor fields that belong to H 3 (U,^) for any relatively compact 
U C M. 

Definition A.l. Let {Uj} be a locally finite covering of M by relatively 
compact open sets Uj . We define (M, 7) as the space of tensor fields 
/ E fff oc (M,7) such that 

II / II fl-;, : = sup || / \\h-(u j , 1 )< oo- 
j 

H^(M, 7) is a Banach space with the above norm. H^(M, 7) will have 
the usual embedding and multiplication properties of Sobolev spaces pro- 
vided all of the -ff s (t/j,7) have them. This will be the case when (M, 7) has 
injectivity radius bounded from below away from zero. In this case, it also 
holds that the Sobolev constants relative to the subsets Ui are uniformly 
bounded, and 7^. is uniformly equivalent to the Euclidean metric. 



Appendix B. Derivation of the reduced system. 

In this appendix, it is shown how equations (3.29) are obtained from (3.14) 
after making a gauge choice. This has been done first by Friedrich in [21], 
which the reader is referred to for more details. 

Assume that a solution to the Einstein-Euler-Entropy system in the 
frame formalism, equations (3.14), is given. 

In fluid source gauge, the vanishing of T a corresponds to equation 
(3.29a), while (3.7), (3.17) and d a MS = give (3.29b). 

It follows from our definitions that 

P > Q 



p- 



dp 
dg 2 



u a u x V \up — upu x S7 \u a 



dv 2 1 /„ r dv 2 \dp 1 9 dg-\ ( _ ,, _ ,, \ 

"5 r - 1 + "2 -7T hf + — ^ V 7T UaV^Sp - UpV^Sa 

os n+ oV v or J os o + v dsl\ / 



ds p + g 
1 (dg I dp 



a o "5~ SaU^V^Ufi - SpU^VnUp 



p + g\ds v 2 ds 

Combining the quantities in the system (3.14) with the above expression, we 
obtain that in fluid source gauge equations (3.29c) and (3.29d) correspond to 

- 7x^v [0 (? a ] = 



P + g 



and 



30,3 p + Q 
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respectively. 

Next, notice that equations (3.29g), (3.29h) and (3.29j) correspond to 
(2.5), (2.15) and (2.8) when written in fluid source gauge, whereas (3.29i) is 
the same as C u \7 a s = 0, which is implied by equations (3.14) (C is the Lie 
derivative) . 

Finally, we have the decomposition 

V^TpQ = wu a + q a 

where w is given by the left-hand side of (2.5). From (3.10) and (3.13), we 
obtain 

~Jtq a = n^7T, T TT a x F aTX TT^ - n^F^uW, 

and 

With the help of these expressions and (3.17) and defining V as in (3.28), it 
follows that 

TT ( ^TT p) u F tlTV u T - ^TT^ir^ir^ir^F^xuZ = u^V^Eap + B^V a «"+ 
E^V^ + V^B u{a e^ - 2u x V x u^ {a B 0> - 3n {a ^ M u x E fj)x 
- 2^ x V x u [a E 0) ^ + 7r Q ^ A V A u^„ + 2TT^Tr^V x u l/ E a p 

which, in fluid source gauge, gives (3.29e). The equation (3.29f) is obtained 
by a similar argument after computing 
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